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Abstract. We study the use of a pair of qubits as a decoherence probe of a non-trivial 
environment. This dual-probe configuration is modelled by three two-level-systems 
which are coupled in a chain in which the middle system represents an environmental 
two- level-system (TLS). This TLS resides within the environment of the qubits and 
therefore its coupling to perturbing fluctuations (i.e. its decoherence) is assumed much 
stronger than the decoherence acting on the probe qubits. We study the evolution 
of such a tripartite system including the appearance of a decoherence-free state (dark 
state) and non-Markovian behaviour. We find that all parameters of this TLS can be 
obtained from measurements of one of the probe qubits. Furthermore we show the 
advantages of two qubits in probing environments and the new dynamics imposed by 
a TLS which couples to two qubits at once. 
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1. Introduction 

The loss of coherence (decoherence) of quantum bits (qubits) due to environmental 
perturbations is an important obstacle on the way to large scale quantum electronics 
and quantum computation. Such perturbations at the same time contain information 
about the surrounding environment which generates them. The idea of using qubits as 
probes of their environment has recently generated interest [H [21 [3l HI [5] as an alternative 
application of qubit technology where the effects of decoherence are used, rather than 
suppressed. 

In general, when an environment acts on a qubit as a weakly coupled, fluctuating 
bath, the environmental effects can be simply expressed as a relaxation and an 
excitation rate as well as a pure dephasing rate. The decoherence process becomes 
much more complex when a qubit couples to any component of an environment 
strongly enough such that quantum mechanical levels within the environment need 
to be taken into account. In many systems, such partially coherent "pockets" are 
observed in the environment. Examples include two- level fluctuators in superconducting 
devices [HI El El El [IHl IH], impurity spins in semiconductors [121 [131 [TH [TSl HE] and 
single-molecule magnets p!7l [T8l 119] such as ^Fe and Ferritin. Valuable information 
about the quantum mechanical nature of the environment may be obtained by using 
qubits as a direct environmental probe, with potentially high sensitivity and high spatial 
resolution 

The concept of qubit probes has already been used to realise a nano-magnetometer 
using Nitrogen- Vacancy (NV) centres in diamond attached to the end of atomic force 
microscope cantilevers [201 [211 [221 [231 [21] ■ In this case the Zeeman splitting of electronic 
levels within the NV centre is optically probed to provide a direct measure of the local 
magnetic field to nanometre resolution [21]. Additional information which one may 
obtain from the decoherence processes has also been studied [251 [26] in this context of 
classically fluctuating fields. 

We investigate a more general configuration in which a decoherence probe couples 
transversally to a localized pocket of coherence and can therefore exchange energy with, 
and extract information from, its immediate environment. We model such a pocket as 
an environmental two- level-system (TLS), for example a charge- or spin- impurity, which 
is then in turn coupled to its surrounding environment. We will show a full solution of 
the system dynamics in the relevant parameter regimes and demonstrate how one could 
use this information to identify and characterize an environmental TLS. 

The behaviour of one qubit interacting with such a partially coherent component 
of the environment is well understood [27] . In such systems ambiguity often arises 
between oscillatory evolution induced by partially coherent defects and oscillatory 
evolution which stems from the single qubit Hamiltonian itself. Two non-interacting 
qubits which couple to the same environment (figure [T]A) enable the observation of any 
environmentally induced correlations between these qubits, removing the ambiguity. 
This indirect interaction is a good indicator for the existence of coherent regions in the 
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closer environment. 

Solving the equations of motion for such a model numerically results in a variety 
of complex behaviours [2HI EHl [301 El] in different regimes of parameter space. In this 
paper we focus on solving the system analytically in key regimes and then use these 
analytical solutions to understand the more general behaviour numerically. 

Our model of two qubits coupled to a common environmental TLS has many 
similarities to the well studied problem of two qubits coupled to a quantum harmonic 
oscillator. This classic Tavis-Cummings model has been widely investigated for its use 
in quantum information processing and studying atom-photon interactions. Interpreted 
in this context our results complement a variety of physical effects which appear in 
the Tavis-Cummings model, including dispersive coupling [32], |33l IMl ESI ESI EZ], ultra- 
strong coupling [SHI EHl SO] and entanglement birth and death [HI |121 SSI Si]. It 
is natural that similar effects appear in both systems due to their formal equivalence 
within the single-excitation subspace. There are however important physical differences, 
as an environmental TLS is typically more localised in space than a quantum harmonic 
oscillator and the symmetry of the TLS allows a more general coupling to its 
environment. 

After introducing the theoretical model in section |2] we present a general analysis 
of the system in section [3] In the weak decoherence regime we find an oscillation of 
energy between the qubits, i.e. an environmentally mediated coupling. We find a 
decoherence-free state (dark state), which leads to formation of stray entanglement 
between the qubits in the decoherence process independent of the decoherence strength 
of the TLS. In our analytical solutions we find a clear threshold between oscillations 
and decay. In order to clarify what constitutes an environmental pocket of coherence 
it is shown that the same threshold divides Markovian from non-Markovian system 
dynamics. We define an effective decay rate of the qubit dynamics which turns out to 
have a linear dependence on the decoherence rates of the TLS in the weak decoherence 
regime and a roughly inverse dependence in the strong decoherence regime. In section |4] 
we interpret our results in the context of dual-probe microscopy and we find that in the 
weak decoherence regime an environmental TLS can be fully characterized and located 
in a substrate. In the strong decoherence regime the TLS can only be located. Section 
[5] puts the theoretical model in context of present experimental qubit realizations. 

2. Model and methods 

In order to study how a qubit probe pair interacts with an environmental TLS we 
construct a simplified model. Consider an experimental setup in which the two qubits 
are attached to an atomic force microscope such that they can be positioned precisely 
(with a fixed distance between them) on top of a scanned substrate which contains the 
TLS (figure [i]A.). At several positions the population of the excited state of the qubit is 
measured as a function of time. In each position of the cantilever the coupling strengths 
between the qubits and the TLS, gi and g2, vary due to their relative position. We will 
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Figure 1. A Experimental setup in which two qubits are attached to an atomic force 
microscope and moved along their connecting line y with a fixed distance between 
them. The TLS is located in the substrate underneath the qubits. B Illustration 
of the model system investigated in this paper: Two qubits with level splitting LOq 
coupled with individual coupling strengths g\ , 52 to an environmental TLS detuned by 
2(5. The TLS in turn is coupled to a bath. C Energies of the Hamiltonian eigenstates 
and the corresponding decohering transitions. Dephasing transitions (oc F^) transfer 
population between states of similar energies, relaxation transitions (oc Fi) transfer 
population to subspaces of lower excitation number. Crosses indicate rates which 
disappear for g\ — gi- The eigenstates are given in appendix Appendix C 



give a detailed model of this variation in Sec. |4j The full system Hamiltonian can be 
written as: 

Hsys = u}q(x'^ +u}ga'^ +[Uq + d)a^ + gi[a^ + <y^ ) 

I „ l„Q2rLS I 02TLS\ /i \ 

+ 92{(r^(T^ +(^y (^y > (1) 

where a^yCTy and are the respective Pauli operator which act on qubit 1 (Ql), qubit 
2 (Q2) or the TLS (TLS). The first two terms describe the two qubits which have the 
same level splitting of 2ug while the third term describes the TLS with a level splitting 
of 2{uq + 6). Here, 6 is the relative detuning between qubits and TLS. The last two terms 
in eq. ([T]) are transversal coupling terms between each of the qubits and the TLS with 
the respective coupling strengths gi and g2 and we introduce g = \/ g{^ g\ for later 
simplicity. We focus on the specific case of transversal coupling as we are particularly 
interested in direct energy exchange between the qubits and the TLS. Throughout this 
discussion, we use the terminology "qubit" and "TLS" to differentiate between the 
fabricated and controllable two-state-probes and the environmental two-level-system of 
interest. For ease of notation we will use the convention H = ks = 1- 
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As the system Hamiltonian Hgys is block-diagonal, the coherent evolution is limited 
to the subspace states with equal excitation number. For the time evolution, we choose 
the state in which qubit 1 is in its excited state and the other two subsystems are in 
their ground state \Q1,Q2,TLS) = \ tii) as the system's initial state. This is a state 
with a single excitation, and therefore we can neglect the subspaces of higher excitation 
numbers in the following calculations. 

A key advantage of probe qubits attached to a cantilever is that they can be 
calibrated while lifted away from the sample. This allows the intrinsic decoherence of 
the probes themselves to be accurately characterized. Once the probes are brought close 
to the sample, this intrinsic decoherence defines a limit in the sensitivity for detecting 
features within the sample. For our purposes, an ideal qubit probe is one with a very 
long intrinsic decoherence time, as this provides a large dynamic range for sensing. We 
assume the intrinsic contribution to decoherence to be small compared to the dynamics 
induced by the environmental TLS. Then we can ignore this contribution in what follows, 
i.e. assume that only the TLS is coupled to a fluctuating bath. For its coupling to the 
environment we take the operator 

Hirrt = sB = [v^al'^'^ + (2) 

where v_[_ and v\\ are the transversal and longitudinal coupling strengths respectively 
and B is an operator acting on the bath. We assume a low temperature bath Uq ^ T. 

In most qubit architectures, the qubit 's level splitting is significantly larger than the 
other energy scales in the problem. Typically this is a requirement to obtain coherence 
over long time scales as well as adequate control over the quantum system. We therefore 
assume throughout this discussion that ujq ^ 6,gi,g2,T. Under this assumption we 
can neglect subspaces with more than one excitation. A large Uq guarantees a clear 
separation of these subspaces in energy while a low temperature bath guarantees the 
absence of spontaneous excitations from the bath. In the limit of large Uq one can 
often make an additional secular approximation which we discuss in detail later on. 
Breaking the assumption that the qubit 's level splitting is the largest energy leads to 
the ultrastrong coupling regime, which is studied elsewhere |45] . 

We model the time evolution of the system's reduced density matrix p using the 
Bloch-Redfield equations [16]' |17]. Using the eigenvectors |1) to |8) (given in appendix 
Appendix C) of Hgys as basis states, the Bloch-Redfield equations read element-wise: 



Pnm ^^nmPnrn ~l~ ^ ^ T^nmn' m' Pn' m' (3) 

n'm' 

with the Redfield tensor: 

T^nmn'm' ■ ^m'mnn' ^nn'm'm ( An/cfen"^mm' ~l~ -^kmm' k^nn'") 

k 

^nmn'm' • ^ nmS n' m' '^C (^LO LOm'n') 
-^nmn'm' ■ ^nm^n'm' (^OJ UJ^i'm'') ■ 
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Here pnm = {n \ p\m) denotes the density matrix element at position n, m and Unm '■ = 
Un — ujm is the energy difference of the Hamiltonian eigenstates \n) and \m) of the 
system. The system operator which couples to the environment s, is defined in Hint- In 
this approach the environment is solely characterized through its spectral function 

/oo 
dT e'"" {B{t)B{0)) (5) 
'OO 

where the bath operator B is taken in the interaction picture defined with respect to 
Hint- This spectral function, eq. |5| is assumed to change slowly such that it does not 
change on the small scale of 6 and g (but only on the much larger scale of 2ujq). 

The transversal coupling to the low temperature bath leads to unidirectional 
population transfers to states with lower excitation numbers, i.e. relaxation. These 
transfers appear in the Bloch-Redfield equations as linear dependencies of the time- 
derivatives of certain diagonal elements of the system's density matrix on other diagonal 
elements, each with a coefficient. These coefficients (i.e. relaxation transition rates) are 
all proportional to vj_C{2ujq)- For later use we define a general relaxation rate due to 
coupling of the TLS to the environment: 

Ti := vlCi2uq) (6) 

The longitudinal bath coupling leads to two processes: First a loss of phase 
coherence between the states of the system, i.e. dephasing, which is mathematically 
represented by the decay of off-diagonal elements in the density matrix. Second a 
mutual population transfer between certain eigenstates |18] with the same excitation 
number: |2), |4) and |5), |7) (figure [l]C). The corresponding decay rates and transition 
rates are all similarly proportional to: 

:= v^^CiO) (7) 

An energy diagram of the eigenstates is depicted in figure [Tp and all transitions are 
shown as arrows. 

When the TLS is decoupled from the qubits {gi = g2 = 0) then Fi is the decay 
rate of the population of its excited state and 2F^ is the additional decay rate of its two 
off-diagonal elements, i.e. its relaxation and pure dephasing rate respectively. 

Comparing the resulting Bloch-Redfield equations with an approach assuming 
Lindblad equations [HI |50] with a phenomenological relaxation rate and dephasing rate 
on the TLS we find that the two sets of differential equations are equivalent when the 
following two conditions are met. First the spectral function should not change on 
the scale of 6 and g- As the second condition one of the following three requirements 
has to be fulfilled: Either (i) the full secular approximation (explained in the next 
paragraph) is applied to both the Lindblad and Bloch-Redfield equations or (ii) we take 
only longitudinal TLS-bath coupling, i.e., f ^ = or (iii) we assume only transversal 
TLS-bath coupling, fy =0, and choose an initial state which is confined to the single 
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excitation subspace. In case of equivalence the two phenomenological rates in the 
Lindblad equations can be identified as our definitions Fi and F^. 

The full secular approximation neglects all dependencies between different elements 
of the system's density matrix if at least one of them is an off-diagonal element. The 
necessary and sufficient condition for this approximation is that the system's level 
splittings and their differences are large compared to the decoherence rates. Physically 
this means assuming cUg ^ (yf ^ Fi, F^^ and ojq ^ \5\ i.e. the TLS is somewhat coherent. 

In the following sections, analytical solutions to the Bloch-Redfield equations in 
different regimes are discussed. For 5 = these solutions are given in appendix Appendix 



[a} The first solution we show (appendix Appendix A. 2 ) is obtained using the full secular 



approximation and is valid for what we call the weak decoherence regime, when the 
resulting decoherence rates are smaller than the coupling strength g between the qubits 
and the TLS. We obtain two further analytical solutions for purely transversal (appendix 



Appendix A. 3) i.e. v\\ = 0, and purely longitudinal (appendix Appendix A. 4) i.e. v± = 0, 
TLS-environment coupling. The combination of our particular initial state, the large 
Uq limit and f_Lf|| = allows us to solve the master equation without the secular 
approximation. That means no assumption about the relative sizes of g and Fi,F2 has 
to be made. These last two solutions are therefore also valid for strong decoherence, 
i.e., when the decoherence rates are bigger than the coupling strength g. 

3. Dynamics 

In this section we present the analytical solutions to the Bloch-Redfield equations for 
our system. We start the section by summarizing the results for a single qubit coupled 



to a TLS for later comparison. Following this, in section 3.2 we study the behaviour of 
two qubits coupled to a TLS in detail. 

3.1. A single qubit coupled to a TLS 

Before we study the more complicated case of two qubits, the dynamics of a single 
qubit coupled to an environmental TLS (see figure [2]A.) provides a clear overview of the 
relevant physics. This special case can be obtained from all solutions by setting g2 = 
and tracing out the second qubit. Performing this on the Hamiltonian, eq. Q, yields 



(for corresponding eigenstates see appendix Appendix C ) 



H,q = uqaf + {uq + + g,{afal^' + a?VJ^^) (8) 

A comparison with the results in this section will later allow us to distinguish 
phenomena which depend purely on the existence of two qubits and those which are 
due to qubit-TLS coupling in general. This section strongly depends on previous work 
[27] . which is reproduced in our notation. As observable we consider the expectation 
value {crf^) (which is proportional to the qubit 's energy). Equivalently, one could use 
the probability to find the qubit in the excited state, Pexc = | ((c?^) + l)- 
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Figure 2. A Illustration of the simplified model system of one qubit with level splitting 
ujg coupled with coupling strength gi to an environmental TLS detuned by 6. The TLS 
in turn is coupled to a bath. This system follows from figure [T]A. by setting 172 — 0. 
B Expectation values for the case of a single qubit coupled to a TLS as a function 
of time for Ti = — O.lgi; g2 = 0; S = C Real part of the one-sided Fourier 
transform of (cr^^)(i) as a function of detuning S with Fi = F;^ = O.I51; 52 = 0. This 
is a numerical solution of either the Bloch-Rcdfield equations, where Fi and F;^ are the 
definitions given by cq. [6] and [7] or a numerical solution of the Lindblad equations with 
phenomenological rates. Analytically we find the angular frequency of the oscillation: 
(jJosc = 2-\/ (5^ + igf which is drawn as a dashed line. For details of the calculation, see 
section 13.11 



Assuming weak decoherence (i.e. we take the full secular approximation in the 
Bloch-Redfield equations) and no detuning 6 = one finds the expectation values as a 
function of time: 

{af){t) = -1 + e-"^' + e^-^-^^^' cos [Ag,t] (9) 
= -1 + e""^* - e("'^~^'')*cos [Ag^ t] (10) 

The population oscillates between the qubit and the TLS (cf. Fig |2]B) with the 
oscillation frequency proportional to their transversal coupling strength. The oscillations 
decay on the time scale corresponding to the decoherence rates of the TLS. 
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Equivalently we can consider the evolution in Fourier space, where the frequency 
and the decay rate are equal to the position and width respectively of the corresponding 



frequency peak (for details see appendix Appendix A.5 ). In figure[2p the real part of the 
one-sided Fourier transform of {o'f^){t) is plotted as a function of detuning between the 
qubit and the TLS. The peak which starts at frequency Agi diminishes with increasing 
detuning, indicating a shift from oscillatory behaviour to pure exponential decay. The 
analytical solution for {o'f^)it) with 6^0 contains complicated coefficients [27], but the 
frequency of the oscillation is simply Uosc = 2a/5^ + 4gf. This corresponds to the level 
splitting between the hybridized states and [4)^^^ (given in appendix Appendix C) 



and is plotted as a dashed line in figure [2p. The oscillatory behaviour is described by 
this one frequency, which corresponds to the standard "generalized Rabi frequency" [51] 
from quantum optics. For large detuning 6 the expectation value {o'f^){t) is dominated 
by one purely decaying term. A Taylor expansion for small g/6 on the corresponding 
decay rate shows that the rate vanishes with increasing 6/g as (Fi + 4:T^)g^/S^ — )• 0. 

So far only weak decoherence on the TLS was considered. We also want to consider 
the limit where the decoherence is stronger than the qubit-TLS coupling. Simplifying the 
equations to purely transversal bath coupling (fy = ^ F^ = 0) one finds analytical 
solutions for the system dynamics without the use of the secular approximation and 
therefore valid for stronger decoherence: 



(af)(t) = -l + 



-64g? tri 

— —e~'~ 



2e- — 



-^2gl + VI 



cosh 





H ^ sinh 


tjji 









:ir 



where /i := ^JT\ — GAgf. From this expression we see that as the decoherence rate 
Fi increases relative to the qubit-TLS coupling strength gi, the dynamics changes from 
oscillations to pure exponential decay. This becomes obvious by rewriting the hyperbolic 
cosine as: 



cosh{-\/Tl-6Aglt 



cos I -^/6Agj~Tjt 



1 



(e+- + e--) 



for 8gi > Fi 
for 8g, < Fi 



(12) 



and similarly for the hyperbolic sine functions. Therefore, we can identify the threshold 
between oscillations and decay in our approximations as precisely Fi = 8gi. 



3.2. Two qubits coupled to a TLS 

Having reviewed the behaviour of a single qubit coupling to an environmental TLS, we 
now consider the behaviour of a dual-probe configuration. Such a system is of particular 
interest when there is no direct coupling between the qubits. This situation allows us to 
probe what we call coherent pockets of the environment. When such a pocket is present 
in the environment, probing simultaneously with two qubits shows qualitatively different 
behaviour to the standard weakly coupled, Markovian environment, which would affect 
each qubit independently. 
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Figure 3. A Expectation values of botli qubits and tlie TLS. Tlie excitation is shifted 
from one qubit via the TLS to the other qubit and back. The parameters are chosen 
as: Fi = = O.lg; gi — g2; 6 = B Real part of the one-sided Fourier transform 
of {(j'^^){t) as a function of detuning S. Analytically we find the three oscillation 
frequencies —S + \/ 6^ + Ag^ , S + \/ + Ag^ , 2y^S^ + Ag"^ which are given by the 
dashed lines in the plot. The parameters are: Fi = Fj^ = O.lg; 5i = <?2 C Time 
evolution of the expectation values for strong detuning 6 = 5g; Fi = F^ = O.lg; gi = 
52- Even though the TLS is only minimally excited, the effective coupling mediated 
by it still leads to coherent exchange of energy between the two qubits. 



3.2.1. Mediated coupling between the qubits in the weak decoherence regime We will 
first show how the coupling between qubits and TLS will mediate an effective interaction 
between the two qubits themselves. For simplicity we initially consider 5 = 0, i.e. both 
qubits are resonant with the TLS. In this case, and for the initial state chosen in section 
|2] the energy from qubit 1 coherently oscillates between the two qubits via excitation 
of the TLS (figure [3]A.). In contrast to the simpler case presented in section 3.1, the 
oscillations now show two distinct frequencies, namely Ag and 2g. The smaller of the 
two frequencies corresponds to the oscillation of energy between the two qubits. The 



full analytical expression for (cr^^) can be found in appendix Appendix A 



Again the change in the system dynamics due to detuning 5 7^ can be understood 
best by regarding the Fourier transform of the expectation value {crf^){t). This yields 
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a peak for each term located at the corresponding frequency whose half-width-at-half- 
maximum (HWHM) equals the decay rate of the corresponding oscillatory component. 
Figure |3]B shows the result of a numerical solution of the Bloch-Redfield equations with 
dashed lines indicating the analytical expressions for the frequency shifts due to the 
detuning. With increasing detuning 6, the peak at frequency 2g splits into two different 
frequency peaks. The amplitude of the two high frequency contributions diminishes 
with stronger detuning 6, while the amplitude of the lower frequency peak increases. 
This means for stronger detuning S > g, the energy oscillates between the qubits mainly 
at the lower frequency ujiow = l*^! ~ — 4(7^. For sufficiently strong detuning, the TLS 
is largely unpopulated during this process (figure |3]C). 

This kind of off- resonant interaction with the TLS leads to an effective transversal 
coupling between the qubits. This is the usual dispersive coupling term [^21 E21 E31 EH 
[35] . in this case due to virtual excitation of the TLS. Performing a Taylor expansion 
for 5f/(5 <^ 1 on both the lower oscillation frequency uiow and on the decay rate of this 
oscillating term in the weak decoherence solution with detuning yields uiow ~ 2gy\6\ and 
liow ~ ^'^(ri + 12r^)/65^. The lower frequency can here be interpreted as the effective 
coupling strength between the qubits uiow = deff- This effective coupling approaches 
zero slower than the decay rate 'yiow as the magnitude of the detuning increases. The 
strongly detuned TLS therefore mediates an effective transversal coupling between the 
qubits with a weakened influence of the TLS' decoherence rates. However, ultimately 
the effective coupling strength (i.e. the frequency of the oscillation) approaches zero for 
S/g — )■ oo. 

Here we see a fundamentally different behaviour as compared to a single qubit 
coupled to a TLS. The oscillations do not change to a pure decay for strong (but not 
yet inflnite) detuning 6. Additionally the frequency of these oscillations approaches zero 
much slower (oc g/6) than the decay rate of the single qubit (oc g'^/6'^). 

This result has important implications for future experimental designs involving 
several qubits in a closely conflned space. There the occurrence of an environmental 
TLS, which couples to two qubits at once might have a realistic probability, especially 
in solid-state qubits. In that case the qubits are affected by the TLS over a wide range 
of detuning, causing effective coupling between the qubits. 



3.2.2. Formation of stray entanglement As we can see in Fig 3 A, the steady state of 
both qubits is not their respective ground state. Rather, they decay into a state with 
a flnite probability of flnding them excited. This behaviour can be attributed to the 
existence of a so called dark state in our system. The state |3) = y ~ ^ liti) 
(appendix [Appendix C ) is an entangled state of both qubits with the TLS in its ground 
state. The amplitudes of the two states (with the respective qubit excited) have a relative 
complex phase of vr in the time evolution which leads to a cancellation of the qubits' 
influence on the TLS. In our system, this state is thus not influenced by decoherence 
and the system will remain in it for a long time (i.e. for the intrinsic decoherence time of 
the qubits). This is simply a manifestation of the physics of super- and sub-radiance [53] 
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due to the interfering pathways from the qubits to the TLS. Since our chosen initial 
state is a statistical mixture including the eigenstate |3), the steady state of the system 
will still include this fraction of the dark state. The entanglement of the two qubits in 
the steady state depends on the interplay of two things: the "concurrence" [51] of the 
dark state by itself which is given by C = 2gig2/g'^ (i.e. a Bell state for gi = g2) and the 
fraction of the dark state in the mixed steady state. Taking both into consideration we 
find the maximal "entanglement of formation" [54j of the final state a.s E = 53% which 
is reached for gi = g2/ VS- 



3.2.3. Threshold between weak and strong decoherence When the decoherence rates of 
the TLS are stronger than the qubits-TLS coupling Fi, > g the secular approximation 
(section [2]) can no longer be fully applied. With increasing decoherence rates the 
dynamics of the three subsystems changes from an oscillating behaviour to a pure 
decay. This behaviour is analogous to a single qubit coupled to a strongly decoherent 



TLS. In section 3.1 we saw that in this case the crossover was defined by the point 
Fi = 8gi. For two qubits, the crossover between the weak and strong decoherence 
regimes is investigated numerically. For oscillations to occur between the qubits and 
the environmental TLS there needs to be an instant in time in which the population of 
the TLS is larger than both qubits combined. We therefore define the maximum value 
in the evolution: 

M = max {{aj^'m - [(af )(t) + (af )(t)] } (13) 

as a measure of the strength of the oscillation. In the regime of strong decoherence the 
energy of the qubits decays via the TLS to the environment and our defined measure 
is always zero. In the weak decoherence (i.e., oscillating) regime the energy leaves 
the qubits and then partially returns via coherent oscillations from the TLS. This 
gives a positive value for the defined measure. Figure |4] is a logarithmic plot of this 
measure as a function of the two decoherence rates Fi and F;^. There is a sudden 
drop in the oscillation strength to negligible values marking a clear threshold between 
the oscillating (weak decoherence) and the decaying (strong decoherence) regime. The 
oscillating regime (dark area) also marks precisely the parameter regime in which the 
full secular approximation is valid. 

For purely transversal (respectively purely longitudinal) TLS-bath coupling = 
^ F^ = (respectively = ^ Fi = 0) an analytical solution can be found. 



Analogous to eq. (11) and (12) we find the analytical threshold between strong and 
weak decoherence regime at the two points: 

Fi =8^,F^ = 

and (14) 
F^ = 4^?,Fi=0 

This corresponds to the point where the threshold in figure |4] crosses the two axes. The 
factor of two between Fi and F^ stems directly from the definition of the rates (eq. 
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Figure 4. Logj^Q of the oscillation strength (eq. 13 1 as a function of the decoherence 
rates Fi and F^. The threshold between oscillations and decay can be seen as a drop in 
the oscillation strength by four orders of magnitude. We solved the full Bloch-Redfield 
equations numerically (without the secular approximation) and for simplicity we set 
gi — 52 and (5 = 0. The level splitting was Uq — lOOOg. 
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TLS 



Reservoir 
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Figure 5. Two possibilities of placing the system/environment boundary: The left 
case was used to set up and solve the Bloch-Redfield equations for which the system 
dynamics is by default assumed to be Markovian. For the configuration on the right 
side, where the TLS is seen as a part of the environment, we investigate Markovianity 
of the system of two qubits, see text. 



|6] and [7]) as they appear in the master equations i.e. the off-diagonal elements of the 
uncoupled TLS-density matrix decay with a rate 'iT^p). 

3.2.4- Markovianity The coherent coupling to environmental states usually leads to 
non-Markovian lini |55] dynamics in the system (excluding the environmental states). 
Using our model, we can choose where to draw the system/environment boundary (see 
figure [5]) and therefore explore this behaviour in a systematic fashion. Regarding the 
TLS as part of the system, Markovian dynamics is assumed by default as this is a 
necessary condition to apply the Bloch-Redfield equations. However, tracing out the 
TLS we investigate Markovianity of the two qubit system. 

Several measures exist to theoretically quantify non- Markovianity |5 6 1 [571 158] in 
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qubit systems. For our purposes we use the connection between the degree of non- 
Markovianity and an increase in the time evolution of the distinguishabihty, calculated 
with the trace-distance D{t) = Tr | A/[pi(i) — P2(i^)]^| of two different initial states. 
By tracing out the environmental TLS in the analytical solution for purely transversal 
TLS-bath coupling (thus reducing the density matrix to the system of the two qubits) 
one can find the trace-distance between the two time evolutions of the initial states 
\Q1,Q2) = and \Q1,Q2) = |it))- We find that the trace-distance only increases in 
the weak decoherence (oscillating) regime. This means that it is the same analytically 
exact threshold between oscillating behaviour and decaying behaviour (compare eq. 



(14)) that divides non-Markovian from Markovian behaviour. The chronologically first 



and strongest increase AD^ in the trace-distance is 



I exp I ^ = 1 for 8(7 > Fi 

AD,= { v/6V-F?i (15) 



for 8^ < F 



for the particular case of pure relaxation and equal qubit-TLS couplings (F;^ = 0, gi = 
g2 = g/\/2), where AD^ = corresponds to a monotonically decreasing trace-distance, 
i.e. no increase in D{t). Here we see that our experimentally measurable definition of 
the oscillation strength Ai is in this case equivalent to a measure of non-Markovianity. 
We can therefore see the direct link between probing coherence in the environment and 
probing non-Markovianity. 

3.2.5. Effective decay rate In order to characterize the decohering influence of the TLS 
on our probe qubits, we want to introduce a single, effective decay rate, serving as a 
figure of merit towards determining the influence of the TLS' decoherence rates Fi,F^. 
In general, the time evolution of the qubits' expectation values are given by sums of 
exponential functions, where, for the weak decoherence regime, some of the terms will 
be decaying oscillations. In the strong decoherence regime, the effective decay rate is 
introduced by replacing the sum of exponentials f{t) = 6^P(~7i^) t>y a single 

exponential cexp(— 7e//t) with the two conditions: 

/•oo /»oo 

1) / dtf{t)= dtcexp{--fefft) (16) 
Jo Jo 

2) /(0) = c (17) 



which leads to the simple formula: 

7e// 



This defines a single decay rate in the strong decoherence regime. From the two 
analytical solutions which are valid for strong decoherence given in appendix Appendix 
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|A](zero detuning and either purely longitudinal or purely transversal TLS-bath coupling) 
we see that 

8g^ {gl + Agj) 

For purely transversal bath coupling the effective decay rate 'yeff,± is plotted in 



figure |6| Interestingly the effective decay rate eq. (19) and (20) monotonically decreases 
with decreasing ratio g/Ti {g/Tip respectively). The observation of weaker decoherence 
on the qubit with stronger decoherence rate of the TLS is due to "blocking" of dynamics 
of the TLS due to strong decoherence: The exchange of energy between the qubits 
and the TLS is slowed down and thereby also the loss of energy to the environment. 
Regarding decoherence as a measurement process this is analogous to the Zeno-effect [IH] . 

This behaviour is in contrast to the weak decoherence regime, where both purely 
decaying exponents and decaying oscillations occur. We find two ways of defining an 
effective decay rate: either describe the decay of the envelope of the oscillations or the 
effective decay of their average. For details, we refer the reader to Ref.[27] and appendix 
Appendix B Here we use the decay of the average and find the effective decay rate 7e// 



to be linearly dependent on the relaxation rate Fi: 

7e// = ^ri forFi<8(7,F^<4(7 (21) 

In the intermediate regime 2 < Fi/^f < 8 and 1 < ^<fi/g < 4, where the lower bound 
is found empirically, the oscillations are slow on the timescale of the decay. In this 
case the first half oscillation, transmission of energy from the qubit into the TLS, 
dominates the behaviour, and the average decay rate does not reproduce the behaviour 
well. An effective decay rate is not a good description of the dynamics in this regime 
and the apparent discontinuity in figure |6] actually appears as a smooth transition in the 
time evolution of the qubits' expectation value. We have also plotted three numerical 
calculations in figure |6] for different level splittings coq (while 5 is always zero). As stated 
in section [2] the analytical result is obtained in the single excitation subspace which 
requires Ug ^ gi,g2,S. We see excellent agreement between the analytical solution and 
the numerical solution of the full Bloch-Redfield equations for Ug > 1000(7. 



4. Probing a single TLS with two qubits: Parameter extraction 

After studying the system and its dynamics in the previous sections, we now interpret 
the results in the context of decoherence microscopy. In particular we focus on the 
ability to obtain the TLS parameters with a dual probe and compare it with a single 
qubit-probe. For this purpose we only consider the Fourier transform of the evolution 
of the qubits' excited-state-population, as this conveniently represents the parameters 
of interest. 
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Figure 6. Effective decay rate 7e// of the energy of qubit 1 as a function of the TLS 
relaxation rate Fi (both rates in units of g). In the weak decoherence regime they are 



hnearly dependent (eq. 21), for strong decoherence the effective decay rate decreases 
with increasing relaxation rate Fi (eq. 



19). The other parameters of the plot are: 



gi = 52, F^ 







coupling strengths 




Figure 7. Characteristic behaviour of the coupling strengths to the TLS in the 
substrate as a function of the position of the two qubits y (compare figure [T|A.). The 
distance between the peaks is controlled by the distance between the qubits dqq. The 
width of the peaks is controlled by the height h above the sample. Here we chose 
dqq — 3h. The coupling strengths are normalized such that the maximum value is 1. 



Now we will give a more concrete form to the theoretical coupling parameters gi 
and g2 from eqjlj If we assume the coupling strengths depend on the distances di and 
d2 between the qubits and the TLS, as gj oc l/dj and the qubits are moved along their 
connecting line above the TLS in the substrate (see figure [T]) then the coupling strengths 
behave characteristically as a function of the position y (figure [T]). 

4.I. Weak decoherence regime - oscillating behaviour 

From an experimental point of view the parameters 6, Fi, F^ of the environmental TLS 
and even the coupling strengths gi,g2 are in general unknown. We first consider the weak 
decoherence regime when the qubits are close enough to the TLS (so that g ^ Fi,F^). 
Then the obtained information of a measurement of {crf^){t) is equivalent to a horizontal 
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Figure 8. Top: Real part of the one-sided Fourier transform of (cr^^)(t) in the weak 
decoherence regime (eq. |A.4 ) for (5 = 0, Bottom: The table shows how the parameters 
could be obtained from a measurement of the plot above. 



line in figure |3p. The positions of tlie tliree peaks give tlie tliree frequencies in figure |3]B 
i.e. tlie necessary information to obtain the level splitting of the TLS 5 and the qubits- 
TLS coupling strength g uniquely. Measuring g at several positions above the sample 
allows the position of the environmental TLS to be obtained from the local minimum 
of g in figure [7] i.e. a single TLS in the substrate can be located. 

The widths and heights of the peaks provide further parameters although they 
have very complicated dependencies. In the case of resonance (5 = 0) however we 
find three peaks which allow an enormously simplified parameter extraction shown 
in figure |8| Experimentally one such plot provides enough information to obtain all 
system parameters: g from the position of the peaks, Fi and from the HWHM, and 
(having obtained these three parameters) gi and g2 from the heights of the peaks. All 
system parameters can be obtained from one measurement of the time evolution of the 
excited-state-population of one of the probe qubits on resonance with the TLS. To reach 
resonance experimentally the qubits could always be tuned to resonance with the TLS, 
once 5 is obtained as explained in the previous paragraph. 

The major difference between the two qubits and a single qubit is the behaviour for 
detuning to the TLS. While the single qubit is effectively decoupled by detuning, the 
addition of a second qubit maintains an oscillating signal via the TLS-induced effective 
coupling between the qubits. As a result strong detuning and weak qubit-TLS coupling 
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Figure 9. Expectation value of qubit 1 for the case of two qubits coupled to a TLS and 
9i = .92 (top) and a single qubit coupled to a TLS (32 — 0, bottom). Both plots show 
two different cases: large detuning S — 4.8g and weak qubit-TLS coupling g = COSFi. 
For all plots = 0. These two different cases can only be clearly distinguished with 
two qubits. 



show two fundamentally different behaviours and can be distinguished with two qubits 
(figure [9]). The system is also sensitive to TLS over a wider frequency range as the 
TLS-induced coupling decreases more slowly with detuning. 

Furthermore the additional two lower frequencies in figure [3^ which correspond 
to oscillations between the two qubits make it possible to obtain the detuning without 
changing the level splitting of the qubits. 



4-2. Strong decoherence regime - decaying behaviour 

Scanning a substrate for isolated TLS one might find very different decoherence strengths 
for each TLS, some of which might be fluctuating so strongly (or coupled so weakly) 
that no coherent oscillations will occur even when the qubits are directly above it. In 
that case the above technique of parameter extraction is no longer apphcable. However, 
the TLS can still be located (both with a single and two qubits) by monitoring the decay 
rate of the qubit at different locations along the y-axis in figure [l]A.. 

For pure relaxation the position dependency of the effective decay rate is shown 



in figures 10 A. (single qubit probe) and [lOp (two qubit probes). The characteristic 



behaviour provides the position of the TLS. 



5. Experimental realisations 

Although in principle any qubit architecture can be adapted for performing decoherence 
microscopy, in order to study microscopic pockets of coherence, atomic scale qubits with 
long coherence times are ideal. In solid-state, this implies spin donors or defects, such 
as semiconductor donors [Sll EHl US] or colour centres in diamond [201 IMl US] • 

As the NV centre in diamond is an experimentally established and well investigated 
system we discuss the requirements to use these centres in a dual-probe configuration. In 
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Figure 10. A Characteristic behaviour of the effective decay rate of the energy of 
a single qubit ((72 = 0) as a function of the position y. This plot is in the decaying 
regime Fi = lOgi and T^p — . B Characteristic behaviour of the effective decay rate 
of the energies of two qubits as a function of the position y. This plot is in the decaying 
regime Fi = lOgi and F;^ = 

recent experiments the intrinsic decoherence of the NV centre is weak and dominated by 
the dephasing|20| 115] . which sets the sensitivity hmit for probing environmental pockets 
of coherence. Isotopic purification of the diamond lattice [21] is currently investigated 
for quantum computation and sensing purposes and will result in much longer coherence 
times leading to a corresponding increase in sensitivity. 

In order to unambiguously identify coupling via the environment, we need to 
minimise or eliminate coupling between the probes. This can be achieved in two 
configurations. Using the nuclear spin states of the Nitrogen within each NV centre 
as the qubit probes provides a strong electron-nuclear coupling to electron spins in the 
environment whilst minimising the inter qubit (nuclear-nuclear) coupling. For both a 
probe-TLS distance and probe separation of 5 nm, we require a T2 > 20 ms to reach the 
probe-TLS oscillation limit. Using current estimates [2H EI] for the dephasing channels 
due to ^'^C spins in the diamond lattice, this requires a ^^C concentration below 0.03%. 
For this probe separation, the inter-probe coupling is a factor of 1000 times smaller than 
the probe-TLS coupling and therefore provides no extra complication to the analysis. 

A second method of achieving strong probe-TLS coupling is to use a pair of NV 
centres whose crystallographic orientation is such that the natural NV-NV coupling 
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is eliminated due to the angular dependence of the dipolar interaction. Although 
there are more serious fabrication challenges with this configuration, the maximum 
dephasing required (T2 > 20/is) is considerably less due to the strength of the electron- 
electron interaction. Such dephasing times are well within the range currently seen in 
experiments using NV centers P^. 

In either of these configurations, detecting sample impurities with large magnetic 
moments such as single molecule magnets (^Fe, Ferritin) [T71 UHl US] is considerably 
easier, even with currently available intrinsic decoherence times. Depending on the 
background field configuration and qubit operating mode, these impurities will induce 
either dephasing dominated or energy exchange processes. 

6. Conclusion 

In this paper we have investigated the concept of a dual-probe decoherence microscope. 
Using a general model, we have studied the key characteristics of such a system 
analytically and numerically. Mapping out the temporal dynamics of the qubit probes 
provides detailed information about a TLS, the simplest example of a pocket of coherence 
contained within the sample. In addition to the TLS' level splitting and the coupling 
strengths to the probes one can obtain its dephasing and relaxation rate which implies 
the coupling strength to its surround environment and that environment's constitution. 
A dual-probe configuration simplifies the measurement process and increases sensitivity 
to detuned TLS. Furthermore we have shown how the oscillation amplitude of 
environmentally mediated coupling between two probes is largely unaffected by detuning 
and decoherence of the mediating TLS, although the frequency of oscillation still depends 
on detuning. The close relationship between environmentally mediated coupling and 
non-Markovian dynamics makes a dual-probe configuration ideal for both probing an 
environment's potential to induce non-Markovian dynamics in a system as well as 
detecting the spatial extent and interrogate pockets of coherence which sit within a 
more complex environment. 
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Appendix A. Analytical understanding of the expectation values and their 
Fourier transforms 

In this appendix we give the three analytical solutions each for no detuning 6 = (i.e. 
6 <^ g). As mentioned before all three solutions were obtained in the subspace of one 
excitation plus the ground state with the assumption Uq ^ gi,g2- Furthermore a low 
temperature bath Uq ^ T was assumed. First the secular approximation is explained, 
then the three solutions are given. 



Appendix A.l. Matrix form of the master equation and the secular approximation 

The Bloch-Refield equations, eq. |3} can always be rewritten as a matrix multiplication. 
To do so, the elements of the density matrix need to be written in a vector. Here we 
choose the particular order: all diagonal elements first, then all off-diagonal elements. 
To solve the equations the resulting Redfield tensor TZ, which is now a matrix, needs to 
be diagonalized. 

The secular approximation means that for this diagonalisation process one can 
neglect off-diagonal elements in the Redfield tensor TZ when there is a large difference 
between the corresponding diagonal elements. In our chosen order we can regard 
separate blocks in the Redfield-tensor (see below). The diagonal elements in the lower 
right block each have a term —iujk, whose magnitude is given by the energy difference 
of the two system states j and k. If these level splittings are large compared to the 
decoherence rates, then the two blocks linking diagonal and off-diagonal density matrix 
elements can be set to zero, i.e. the upper right block and the lower left block. If the 
differences of the energy differences Ujk — ooim are also large compared to the decoherence 
rates, then all off-diagonal Redfield tensor elements in the lower right block can also be 
set to zero. This is what we call the full secular approximation. The resulting Redfield 
tensor is given by: 
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Mathematically the secular approximation is analogous to the rotating wave 
approximation: Separating the Redfield tensor into a coherent part (i.e. the —iujjk terms) 
and a decoherent part we can define an analogous "interaction picture" for the vector 
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of density matrix elements p: 



d ^ 



dt 



coh4 



(A.l) 
(A.2) 

(A.3) 



This results in a new Redfield tensor TZ in this "interaction picture" where all off- 
diagonal elements are multiplied by a rotating term with the rotation frequency equal 
to the difference of the corresponding diagonal 7?.-elements. The secular approximation 
can then be justified analogously to the rotating wave approximation. 

Appendix A.2. Weak decoherence regime ^ Fi, 

The two qubits couple to the same environmental TLS strongly compared to the 
decoherence of the TLS. The full secular approximation is apphed. 

,4 „4 0^2 „2 ^4 ^ , /I ^2 ^2 ^ , ^4 



9 9 9 9 

i-rm = + ^e-^ - ^e(-^-^)*cos[2, t] + ^e(-^-^)* cos [4, t] 



9 9 9 

,2 ^ , ^2 



9 



= -1 + £ie-^ - %e-^-^-*cos[45t] 



Appendix A.3. Purely transversal TLS-bath coupling vn — 



(A.4) 



F^ = 



The TLS-bath coupling is purely transversal, i.e. there is relaxation only. The secular 
approximation need not be applied here, due to our particular choice of environmental 
coupling operator. Therefore this solution is also valid for strong relaxation. 
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where fi := a/F^ — GAg'^ 



Appendix A. 4- Purely longitudinal TLS-bath coupling v± = ^ Ti = 

The TLS-bath couphng is purely longitudinal i.e. there is dephasing only. Again the 
secular approximation need not be applied here, due to our particular choice of initial 
state and environmental coupling operator. Therefore this solution is also valid for 
strong dephasing. 



^ -2glgl + gl 
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q2 / 4„2p 2p N 
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where /xa := y^T^ - Ag^ 
Appendix A. 5. Fourier transform 

The expressions in all expectation values contain constant and oscillating terms, with 
associated decay rates. Terms of this form are more easily understood in the Fourier 
domain. 

For measured signals of the form e~°*cos[6t] the real part of its one-sided Fourier 
transform yields two Lorentzian peaks at the position of plus and minus the frequency 
h and with a half width at half maximum (HWHM) which equals the decay rate a: 

Fitting such peaks allows us to experimentally obtain the frequency and decay rate in a 
precise and simple way, as displayed in figure [8j Additionally, the close correspondence 
to the parameters in the Fourier domain helps to depict frequencies and decay rates at 
the same time in figures |2p and|3p. 
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Appendix B. Calculation of the effective decay rate of a sum of decaying 
oscillations 

In the strong decoherence regime the oscillations are a sum of several exponentials 
Cj exp(— 7jt). To find one effective decay rate we can simply use eq. 18 This 



procedure is sensible when the different decay rates are not too far (i.e. not orders of 



magnitude) apart. Note that eq. 18 can also be calculated from an integration: 



In the weak decoherence regime on the other hand we have additional oscillations 
for several terms, i.e. an expression of the form 

Cj exp{—'yjt)cos{ujt) (B.2) 

j 

where some Uj might be zero and the cosine function might be replaced by a sine function 
for some terms. In such a case (as for example displayed in figure or [sjA.) one needs 
to decide to either take the average or the envelope of the oscillations. The effective 
decay rate of the average neglects the oscillating terms completely and can therefore 
become zero when there are no purely decaying terms in the expression. On the other 
hand the average is unambiguous while the upper envelope and the lower envelope can 
lead to different effective decay rates. This is the reason why the average was chosen in 
figure [6] for the weak decoherence (oscillating) regime. 

The calculation of the average is performed by neglecting all oscillating terms 



and calculating eq. 18 from the rest. Numerically that is easily done by rewriting 
all oscillations in eq. B.2 as exponentials, which yields an expression of the form: 

Cj exp((-7i + i(^j)t) (B.3) 

j 

Then all terms with a non-zero imaginary part in the exponential rate can be neglected 
and the effective decay rate can be calculated as: 

average: 7e// = ^ ^ ,^ (B.4) 

22k Ckhk 

V ' 

pure decays 

where k sums over all purely decaying terms. 

The calculation of the envelope is performed by setting all oscillations (including 



the algebraic sign) to 1 (upper envelope) or -1 (lower envelope). Afterwards eq. 18 
can be applied to all terms. Numerically that can easily be performed by taking the 
magnitude of the coefficients and real parts of the rates for all oscillating terms: 

upper envelope: 7e// = ^k^k + Ya I^'I 5^ 
lower envelope: 7e// = „ ^k k E; I i\ — (B.6) 
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where k sums over all purely decaying terms and I sums over all oscillating terms. For 
the purely decaying terms it is important not to take the magnitude of the coefficients 
in case negative coefficients Ck < occur. 

In principle all terms with a non-zero imaginary part of the exponential rate are 
oscillations. However, when this imaginary part (which is the angular frequency of the 
oscillation) is smaller than the real part (which is the decay rate) then this term decays 
strongly before the time period of one oscillation, i.e. the term looks like a pure (non- 
exponential) decay. For a numerical criterion whether a term should be categorised as 
oscillating or purely decaying one could therefore measure the imaginary part relative 
to the real part for each individual term. However for simplicity of the criterion we 
categorize all terms with an imaginary part of the exponential rate below 0.1 (where 
5^ = 1) as purely decaying terms in our numerical calculations for figure [6j This is about 
one order of magnitude less than the decay rates plotted in figure |6| 



Appendix C. Hamiltonian eigenstates of the system 



For our system of two qubits coupled to one TLS the unnormalised eigenstates indicated 
in figure [Tp are: 

8) = I nt) 

7) = {-5 + V52 + V)| tn) + 2g2\ nt) + 2(71 1 itt) 

6) = -^7iint) + ^72i;tt) 

5) = (-5 - ^/¥+Ag^)\ tti) + 2^72! tit) + 2g^\ ttt) 
4) = 2(71 1 til) + 2(721 m) + (5 + v/5^+V)l lit) 
3) = -(72|tii)+^7i|iti) 

2) = 2(71 1 til) + 2(721 m) + y^^+4?)| tit) 

1) = I ;;;) 

where t indicates an excited state and I a ground state of the two qubits and the TLS 
in the order \Q1,Q2,TLS). 



(C.l) 



Setting (72 = and tracing out the second qubit yields the system of only one qubit 
coupled to a TLS discussed in section |3.1[ Performing these operations on the above 
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states we find: 



|8)^ 

|6)^ 

|5)-> 

|3)^ 
|2)^ 



8)1, = I tt) 

7)1, = {-5 + a/5^+4^)| U) + ^9i\ it) 
6)1, = -^il tt) 

5)1, = {-6 - ^5^+4^)1 U) + 29i\ it) 
4)1, = 2g,\ U) + ^/P + Ig~')\ tt) 
3)ig = fi-il it) 

2)i, = 2(?i|tt) + (5- V52 + 4^2 
1)1, - I tt) 



(C.2) 



tt) 



Several states become equivalent: 

|3)i, = 5i|i; 



19 



2gi 



|5)i.= 



\7K 
|2)i. 



(C.3) 
(C.4) 



(C.5) 
(C.6) 



|6)i, = -^i|8)i. 

We therefore need (as one should expect) only four states to describe this reduced 
system. 
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